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ABSTRACT 

Context. The identification of the carriers of the diffuse interstellar bands (DIBs) is the longest standing problem in the study of the 
interstellar medium. Here we present recent UV laboratory spectra of various polycyclic aromatic hydrocarbons (PAHs) and explore 
the potential of these molecules as carriers of the DIBs. Whereas, in the near IR range, the PAHs exhibit vibrational bands which 
are not molecule-specific, their electronic transitions occurring in the UV/vis provide characteristic fingerprints. The comparison of 
laboratory spectra calibrated in intensity with high S/N observational data in the UV enables us to establish new constraints on PAH 
abundances. 

Aims. From a detailed comparison of gas-phase and Ne-matrix absorption spectra of anthracene, phenanthrene, pyrene, 2,3- 
benzofluorene, benzo[ghi]perylene, and hexabenzocoronene with new interstellar spectra, we aim to infer the abundance of these 
PAHs in the interstellar medium. 

Methods. We present spectra of PAHs measured at low temperature in the gas phase and in a Ne matrix, and present methods to 
derive absolute absorption cross sections for the matrix and gas phase spectra. We have obtained high signal to noise (S/N > 100) 
absorption spectra toward five lines of sight with reddenings of £ B -v = 1-1-6 mag. The spectra cover the 3050 - 3850 A wavelength 
region where the PAHs studied here show prominent absorption features. 

Results. From the observations, we infer upper limits in the fractional abundances of the PAHs studied here.Upper limits in the col- 
umn densities of anthracene of 0.8 - 2.8 x 10' 2 cirr 2 and of pyrene and 2,3-benzofluorene ranging from 2 - 8 x 10 12 cirr 2 are inferred. 
Upper limits in the column densities of benzo[ghi]perylene are 0.9 - 2.4 x 10 13 and 10 14 cirr 2 for phenanthrene. The measurements 
indicate fractional abundances of anthracene, pyrene, and 2,3-benzofluorene of a few times 10~ 10 . Upper limits in the fractional abun- 
dance of benzo[ghi]perylene of a few times 10~ 9 and of phenanthrene of few times 10~ 8 are inferred. Toward CPD -32° 1734, we 
found near 3584 A an absorption line of OH + , which was discovered in the interstellar medium only very recently. 
Conclusions. The fractional abundances of PAHs inferred here are up to two orders of magnitude lower than estimated total PAH 
abundances in the interstellar medium. This indicates that either neutral PAHs are not abundant in translucent molecular clouds, or 
that a PAH population with a large variety of molecules is present. 
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1. Introduction 

The long-standing mystery concerning the nature of the carri- 
ers of the diffuse interstellar bands (DIBs) has received con- 
siderable and renewed attention in recent years after labora- 
tory spectra of potential candidates have become available. 
Most scientists in the field support the general idea that the 
DIBs arise from electronic transitions of gas-phase molecules. 
Recent examples include the 2?, 2i4i, and 2^ transitions in the 
B 1 B l «- X'Ai system of /-C 3 H 2 (Maier et al. 2010), the (0,0) 
band of the A 2 n„ - X 2 Tl g system of the HC4H + diacetylene 
cation (Krelowski et al. 2010b), or the anthracene and naphtal- 
ene cations (Iglesias-Groth et al. 2008, 2010). Similarly, a recent 
theoretical work established an excellent agreement between the 
Ojj band of the l B x -X l A' electronic system of CH 2 CN with the 
8037 A DIB (Cordiner & Sarre 2007). 

DIBs are generally observed in absorption against the con- 
tinua of reddened supergiants in the Galaxy and in extragalactic 
sources. The appearance of DIBs in translucent molecular clouds 



characterized by a wide range of physical parameters indicates 
that DIBs arise from chemically stable carriers. More than 400 
DIBs have been identified toward a single galactic line of sight 
(Hobbs et al. 2009). The individual absorption profiles and line 
strength ratios of DIBs may vary significantly from one sight- 
line to another (Galazutdinov et al. 2008). While the equivalent 
widths of several DIBs correlate well with the reddenings £b-v 
of the background stars, the majority of the DIB carriers appear 
to be enhanced in the outer regions of translucent clouds (Cami 
et al. 1997; Sonnentrucker et al. 1997). Observations of DIBs in 
the Magellanic Clouds (Welty et al. 2006) have shown that some 
DIBs are significantly weaker than toward Galactic lines of sight 
with similar reddening, suggesting that the abundance of the re- 
sponsible carriers scales with the metallicity of the environment. 
General scaling relations have not been established, however, as 
the equivalent widths of the family of C2-DIBS, for example, are 
as strong as in the Galaxy, for lines of sight with comparable 
reddening (Welty et al. 2006). 
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The vibronic progressions established for several DIBs 
(Duley & Kuzmin 2010) point toward an origin from a torsional 
motion of pendent rings, suggesting that the carriers are 'floppy' 
polycyclic aromatic hydrocarbon (PAH) structures such as 9- 
phenylanthracene, parasubstituted biphenyls, tetracene deriva- 
tives, bianthracene, tolane, etc., rather than compact molecules. 
The presence of PAHs in the interstellar medium is now gen- 
erally established, and there is unequivocal evidence that the 
unidentified near-infrared emission bands arise from PAHs 
(Leger & Puget 1984; Allamandola et al. 1985; Tielens 2008; 
Bauschlicher et al. 2010). Mixtures of PAHs or nano-sized hy- 
drogenated dust particles have been proposed to explain the 2175 
A UV bump seen toward stars reddened by translucent material 
(Joblin et al. 1992; Mennella et al. 1998; Steglich et al. 2010), 
and the presence of carbonaceous dust particles in translucent 
material is well established (Henning & Salama 1998). The first 
detection of the Ceo and C70 fullerenes in the young planetary 
nebula Tc 1 has recently been reported by Cami et al. (2010). 

With the present study, we have undertaken a new effort to 
search for PAHs in translucent molecular clouds. While ion- 
ized PAHs exhibit absorption features in the visible and near 
IR, making them attractive carrier candidates for the DIBs 
(Salama et al. 1999), the electronic transitions of their neu- 
tral counterparts are shifted to the UV. Since neutral PAHs 
can be prepared more easily with sufficient density in the gas 
phase and because observational spectra in the UV can be ob- 
tained with high S/N, we focused our interest on this class 
of molecules. In this paper, we present results on our search 
for anthracene, phenanthrene, pyrene, 2,3-benzofluorene, and 
benzo[ghi]perylene. These molecules were measured in the gas 
phase with our supersonic jet/CRDS apparatus and show promi- 
nent absorption features at wavelengths shorter than 4000 A. 
We have obtained high signal-to-noise ratio (S/N > 100) spec- 
tra toward heavily reddened supergiants in the 3050 - 3850 A 
region using UVES. We developed techniques to calibrate the 
CRDS spectra against matrix and solution spectra and to derive 
absolute absorption cross sections for vibronic bands measured 
in the gas phase. We also searched for hexabenzocoronene, for 
which we only have Ne and Ar matrix spectra. We present a 
method to estimate gas phase transition wavelengths from the 
matrix spectra. These methods are described in detail in Sect. 2, 
which also provides the measured gas phase absorption wave- 
lengths and the estimated absorption cross sections of the PAHs 
studied here. The astronomical spectra, obtained toward early- 
type supergiants with reddenings up to £b-v = 1-6 mag, are 
presented in Sect. 3. A detailed comparison with the laboratory 
results is given in Sect. 4.1, and Sect. 4.2 presents upper limits 
of PAH column densities toward the lines of sight studied here. 
Section 4.3 contains an analysis of various diatomic molecular 
absorption lines which are detected toward the observed stars, 
and which are used for consistency checks with previous re- 
sults. The inferred upper limits and the corresponding maximal 
fractional abundances of anthracene, phenanthrene, pyrene, 2,3- 
benzofluorene, benzo[ghi]perylene, and hexabenzocoronene are 
compared to expectations of PAH abundances in translucent ma- 
terial in Sect. 5.2. 



2. Laboratory Data 

The confirmation of potential DIB carrier molecules requires 
laboratory spectra that are obtained under conditions that match 
as close as possible the conditions in interstellar clouds. Such 
conditions - collision-free environments and low temperatures - 



prevail in the expansion of supersonic jets of rare gases seeded 
with the molecule to be studied. In combination with cavity 
ring-down laser absorption spectroscopy (CRDS) (O'Keefe & 
Deacon 1988), low-temperature gas phase absorption spectra are 
obtained, providing the ultimate data which can be compared 
with astrophysical observations. For several reasons, however, 
supersonic jet CRDS spectra are difficult to obtain. It is for in- 
stance difficult to transfer complex molecules into the gas phase. 
In such cases, another low-temperature technique, matrix isola- 
tion spectroscopy (MIS) (Bondybey et al. 1996), is performed. 
In MIS experiments, the sample molecules are dispersed in a 
transparent matrix of rare gas atoms such as Ne or Ar kept at 
cryogenic temperatures of 6.5 K (Ne-matrix) and 12 K (Ar- 
matrix). Compared to laser spectroscopy in supersonic jets, the 
MIS technique has the advantage that wide ranges of wave- 
lengths are scanned at once and that only small amounts of the 
sample molecules are required in general. The main disadvan- 
tage of MIS experiments is that the molecules are no longer 
in a collision-free environment. Interactions with the surround- 
ing matrix material induce shifts of the transition energies and 
a broadening of the absorption profiles. Spectra obtained from 
MIS experiments are thus of limited use in the astrophysical 
context. On the other hand, CRDS experiments do not gener- 
ally allow to obtain absorption cross sections, as it is difficult to 
measure the density of the relevant molecules in the supersonic 
jets. This is more easily achieved for MIS spectra which can be 
readily compared with solvent spectra, where the concentration 
of the dissolved molecules can be accurately measured. We have 
developed methods to derive absolute absorption cross sections 
for supersonic jet spectra and to determine gas phase absorption 
wavelengths from MIS spectra, and present them in Sects. 2.1 
and 2.2 below. 

From the various PAHs studied in our laboratory, we have 
selected six species revealing prominent absorption features in 
the 3050 - 3850 A region, and for which CRDS spectra are 
available (Rouille et al. 2004, 2007; Staicu et al. 2004, 2008). 
In the order of increasing size, these molecules are anthracene, 
phenanthrene, pyrene, 2,3-benzofluorene, benzo[ghi]perylene, 
and hexa-pen'-hexabenzocoronene (hexabenzocoronene in the 
following). The structures of the six molecules are illustrated 
in Fig. 1. Calibrated solution spectra obtained in the nonpo- 
lar solvent cyclohexane for anthracene, phenanthrene, pyrene, 
and benzo[ghi]perylene (Karcher et al. 1985) and in 1,2,4- 
trichlorobenzene for hexabenzocoronene (Hendel et al. 1986) 
are used as references for the determination of absolute absorp- 
tion cross sections as described in the next section. The solution 
spectrum of 2,3-benzofluorene has been measured in our labo- 
ratory, using cyclohexane as solvent. Spectra of anthracene and 
pyrene have been measured as well and compared with the litera- 
ture data to validate our experimental setup. The matrix isolation 
spectra of the six molecules were recorded in Ne at 6.5 K with an 
apparatus described elsewhere (Steglich et al. 2010). For hexa- 
benzocoronene, the calibration procedure has been restricted to 
a Ne matrix spectrum (Rouille et al. 2009) since a jet spectrum 
in the wavelength range of interest is not available. 

2.1. Determination of absolute absorption cross sections for 
supersonic jet spectra 

In cavity ring-down spectroscopy (CRDS) experiments, the de- 
termination of vibrational absorption cross sections is generally 
not possible, as the determination of the densities in the super- 
sonic jets remains difficult. In contrast, absolute absorption cross 
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Fig. 1. Schematic carbon structures of the PAHs discussed in this 
study, a): anthracene, C14H10; b): phenanthrene, Ci4H 10 ; c): pyrene, 
C 16 H 10 ; d): 2,3-benzofluorene, C 17 H 12 ; e): benzo[ghi]perylene, C22H12; 
f): hexabenzocoronene, C42Hi 8 . 
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sections are routinely obtained from solvent spectra (Karcher 
et al. 1985), where the absorption is expressed in terms of a mo- 
lar extinction coefficient e (in 1 mol 1 cm 1 ). Indeed, assuming 
that absorption is the only significant contribution to the extinc- 
tion, e can be converted into an absorption cross section <x (in 
cm 2 ) according to 



1000 In 10 



N A 



-e= 3.8235 x 10 _21 e 



(1) 



where Na is the Avogadro constant. By integrating cr over a 
wavelength interval that includes a specific absorption band, one 
obtains the integrated cross section of the band. It can be used to 
derive molecular column densities from interstellar spectra using 
the formalism summarized in Sect. 4.2 (in particular Eq. 5). 

Solvents modify the band intensities of the dissolved 
molecules and cause shifts of the transition wavelengths and a 
broadening of the vibronic bands. The effect of the wavelength 
shifts will be addressed in detail in Sect. 2.2. Due to the broad- 
ening, an absorption peak measured in solution may correspond 
to several vibronic bands. However, if solution spectra are suf- 
ficiently well-resolved, they may be used to derive estimates 
of the vibrational absorption cross section for matrix spectra 
(Ehrenfreund et al. 1992) or directly for jet spectra (Kokkin et al. 
2008). The method is based on the reasonable assumption that 
the absorption cross section of a molecule integrated over an en- 
tire electronic transition is independent of the environment of the 
molecule, i.e. whether it is vacuum, a rare gas matrix, or a solu- 
tion. The error that we introduce when we make this assump- 
tion is discussed in the next paragraph. For the calibration of our 
high-resolution gas phase spectra, we have developed a proce- 
dure where we use a matrix spectrum as an intermediate step 
between the solution spectrum and the jet spectrum. We inte- 
grate the complete region corresponding to an electronic transi- 
tion in the solution and matrix spectra and equal them. The same 
operation is then repeated for the single band of interest in the 
matrix and jet spectra. The calibration procedure is exemplified 
in Figs. 2 and 3 for the case of 2,3-benzofluorene. Both figures 
show the measured absorption profiles of 2,3-benzofluorene in 
a Ne matrix, in a solution, and in a jet. The different shifts of 
the band groups seen in Fig. 2 allow us to distinguish between 
the different electronic transitions of the molecule. The calibra- 
tion has been applied exclusively to the first electronic transition 
(Si <— So) of 2,3-benzofluorene. The direct comparison of the 
areas of the higher-energy transitions in the solution and matrix 
spectra is actually complicated by the presence of a continuum 
in the matrix spectrum attributed to the scattering by the PAH 



Fig. 2. Absorption spectra of 2,3-benzofluorene measured in a super- 
sonic jet (red curve), in a Ne matrix at 6.5 K (green curve), and in cy- 
clohexane at room temperature (blue curve). Note that the matrix and 
solution spectra were magnified by a factor of 20. All spectra were mea- 
sured in our laboratory. 



molecules themselves. In Fig. 3, the spectra of the first elec- 
tronic transition are displayed in wave number, and the matrix 
and solution spectra are shifted to the blue to account for envi- 
ronment effects (34 and 491 cm -1 , respectively). To calibrate the 
matrix spectrum against the solution spectrum, the integration 
has been performed from 29 500 to 31 500 cirT 1 (see Fig. 3). 
The area of the gas phase origin band at 29 894 cirT 1 has then 
been equaled to the corresponding band at 3348 A in the matrix 
spectrum while the contribution of the smaller band at 29 872 
cirT 1 has been subtracted. 

In order to estimate the error that we introduce with the 
assumption that the integrated cross section is not affected by 
the environment of the molecule, we have to consider the elec- 
tronic interaction between the solute molecule and the neigh- 
boring solvent molecules. The variations of the band intensi- 
ties caused by this interaction are described by various models 
(Chako 1934; Linder & Abdulnur 1971). In the first approxima- 
tion, the Lorentz-Lorenz force which arises from the polariza- 
tion of the surrounding solvent molecules can be described in 

terms of the Chako factor 9n / (n 2 + 2) , where n is the refrac- 
tive index of the solvent. Multiplying the oscillator strength of a 
band measured in solution by the Chako factor, the correspond- 
ing oscillator strength in the gas phase is obtained. These models 
have been successfully applied to C60 by Smith (1996). For the 
wavelength range and the experiments discussed here, where cy- 
clohexane was used as solvent, the Chako factor is 0.78. In view 
of the fact that the Chako factor is still an approximation and that 
the proper ratio of the oscillator strengths in the gas phase and 
in a solution depends on the exact nature of both the solvent and 
solute molecules, we chose not to apply any correction factor to 
our cross section values. Taking into account the other sources of 
errors including the experimental and data processing uncertain- 
ties, we estimate that the oscillator strengths are accurate within 
±50 %. 

The results of the calibration procedure are summarized in 
Table 1. The wavelength positions and widths are from our 
CRDS experiments, except for hexabenzocoronene, for which 
the wavelength position in the gas phase was extrapolated from 
matrix measurements (cf. Sect. 2.2). Although the band posi- 
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Table 1. Spectral properties (attribution, position, width, integrated cross section, and oscillator strength) of the PAH bands considered for com- 
parison with astrophysical data in this article. 
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Notes. {a) Band positions are given for air under standard conditions (see text). {b) An upper limit of the column density has been calculated 
assuming a detection limit with an equivalent width of 1 mA. Staicu et al. (2004). (<,) Rouille et al. (2010). (c) Rouille et al. (2004). <0 Staicu et al. 
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Fig. 3. S i «— So absorption spectra of 2,3-benzofluorene under the same 
experimental conditions as in Fig. 2. The solution and matrix spectra 



have been shifted to higher energy by 491 cm 
tively. 



and 34 cm , respec- 



tions were originally determined for vacuum, we converted them 
to air because high-resolution astronomical spectra are tradition- 
ally calibrated against Th/Ar lines whose wavelengths are given 
for air. The conversion was carried out according to /t a ; r = A vic /n 
with n being the index of refraction of air under standard con- 
ditions at /l vac , as computed with Eq. Al of Livengood et al. 
(1999). For anthracene and phenanthrene, the results are given 
for the origin of their respective S\ <— S o transitions. For pyrene, 
the spectrum of the S 2 <— So transition is composed of band 
clusters due to the vibronic interaction between the S 2 and S 1 
states. We give in Table 1 the central wavelengths and the to- 
tal FWHMs of the two band groups at 3166 A and 3208 A 
(Rouille et al. 2004). For 2,3-benzofluorene, the origin and a 
strong vibronic band of the Si <— So transition are listed. For 
benzo[ghi]perylene, the spectral region around 3500 A shows 
two broad bands attributed to vibrationally excited bands of the 
S 2 So transition. We also list the oscillator strength / ob- 
tained in the approximation of weak bands (Mulliken 1939) from 
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4eoOTC 

e 2 A 2 



s 



cr(A)dA, 
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where A is the central wavelength of the band and eo, c, m, and 
e are the vacuum permittivity, the speed of light in vacuum, the 
mass, and the charge of the electron, respectively. It is interest- 
ing to note that, in some cases, the absorption cross sections and 
oscillator strengths vary significantly for the different molecules 
studied especially if we consider the transition with the lowest 
energy (Si <— So)- Thus, we found that the oscillator strengths 
of this transition differ by two orders of magnitude if we compare 
the two molecules anthracene and phenanthrene, which are both 
composed of three aromatic cycles (see Figure 1). This shows 
that the change of the electronic structure during the transition, 
which determines the strength of the transition, does not obey 
simple empirical rules. In fact, laboratory experiments seem nec- 
essary to determine this important property. 

In the case of 2,3-benzofluorene, an oscillator strength of 
0.0085 is derived for the origin band of the Si <— So tran- 
sition (see Table 1). The summation over the complete elec- 
tronic transition (i.e. the origin band and the associated vi- 
bronic bands) yields an integrated oscillator strength of 0.03. 
This value is compared in Table 2 with the results of theoretical 
calculations (Staicu et al. 2008) and a laboratory measurement 
(Banisaukas et al. 2004). Theoretical transition energies and os- 
cillator strengths were obtained by applying different quantum 
chemistry models. In Table 2, we refer to results obtained with 
a semiempirical model, ZINDO, an ab initio method, CIS(D)/6- 
31G(d), and a model based on the time-dependent density func- 
tional theory, TD-DFT-B3LYP/TZ. For details on the computa- 
tional methods, see Staicu et al. (2008) and references therein. 
As far as the oscillator strength is concerned, the largest devia- 
tion from experimental data is found for the ZINDO calculation 
although it yields the best result for the transition energy. The 
oscillator strengths of the other two calculations come close to 
the experimental values. The satisfactory agreement of our os- 
cillator strength with ab initio calculations and the result derived 
from Ar matrix data gives us confidence that our method to de- 
termine absolute absorption cross sections leads to reasonable 
results. 

For the other PAHs discussed in this paper, absorption cross 
sections of calculated spectra are available in a database (Malloci 
et al. 2004, 2007). Unfortunately, the theoretical spectra were 
convoluted with a too broad line shape function making a com- 
parison with our experimental data very difficult or even impos- 
sible. 

At this point, a general statement about computed electronic 
transition energies seems appropriate. As can be inferred from 
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Table 2. Comparison of the energies and oscillator strengths for the first 
electronic transition of 2,3-benzofluorene in calculated and experimen- 
tal spectra. 
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" Staicu et al. (2008) 

b Banisaukas et al. (2004) 
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the data of Table 2, the theoretical values may strongly differ 
from the experimental results. The current quantum chemistry 
models and calculation techniques allow one to predict elec- 
tronic transition energies with an accuracy of the order of 0.2 eV 
at best. The accuracy of semiempirical calculations, even if the 
models are well parametrized, suffers from their simplicity. On 
the other hand, the ab initio approach is the most complete. In 
order to be applied to large molecules, however, approximations 
have to be made, affecting the accuracy of the results. Finally, 
calculations based on the density functional theory currently 
miss functionals adapted to PAH molecules. For our example 
molecule, 2,3-benzofluorene, the best result is obtained with the 
semiempirical ZINDO method. However, the deviation from the 
experimental value is with more than 20 A still quite significant. 
In comparison, the extrapolation of matrix data, an experimen- 
tal method which will be described in the next section, yields a 
frequency position as close as 5 A to the gas phase position. In 
conclusion, we can say that the uncertainty of quantum chemical 
calculations is still too large for a meaningful comparison with 
high-precision astronomical observations. 

2.2. Determination of gas phase absorption wavelengths 
from MIS spectra 

PAH absorption wavelengths obtained from MIS experiments or 
from solutions suffer from a broadening and a shift of the ab- 
sorption features, compared to the absorption profiles obtained 
in an environment free of interactions. It is nevertheless possi- 
ble to exploit MIS results to predict band positions in the gas 
phase if they are not available, as in the case of the /3-band of 
hexabenzocoronene (Rouille et al. 2009). For PAH molecules 
in rare gas matrices, the dominating mechanism in the interac- 
tion between the molecules and the rare gas atoms is the dis- 
persion effect. The effect arises from the interaction of polar- 
izable species and is proportional to the product of the mean 
static dipole polarizabilities (London 1937). As the dispersion 
effect depends on the polarizability of the matrix material, a 
given transition shows different energy shifts in different ma- 
trices, and an extrapolation to a matrix with zero polarizability 
allows one to estimate the band position in the gas phase. The 
models of Longuet-Higgins & Pople (1957) and Shalev et al. 
(1991a,b, 1992) have been applied to evaluate the transition en- 
ergy shifts for PAH molecules surrounded by rare gas atoms 
(Biktchantaev et al. 2002; Shalev et al. 1991a,b, 1992). These 
models take into account the arrangement of the rare gas atoms 
around the molecules. As this information cannot be obtained 
by straightforward methods, we assume that the dispersion ef- 
fect does not vary significantly when we consider the different 



geometries encountered for PAH molecules in Ne and Ar matri- 
ces. The transition frequency v gas p h a se of a band is then obtained 

as 

^gas phase = ^Ne matrix + „ 7 (^Ne matrix — ^Ar matrix) > (3) 
K a - 1 

where R a = 4.13 is the ratio of the polarizabilities of Ar and Ne 
(Radzig & Smirnov 1985). 

We apply this method to the complete electronic spectrum 
of hexabenzocoronene that we have recently measured in a Ne 
matrix (Rouille et al. 2009). The full spectrum reveals strong 
absorption bands in the 3000 - 4000 A spectral region. A well 
defined feature of the a-band was found at 4383 ± 5 and 4344 
+ 5 A in Ar and Ne, respectively. The linear extrapolation based 
on the ratio of the rare gas polarizabilities leads to a wavelength 
of 4332 ± 9 A. This feature was measured in a cold molecular 
beam (Kokkin et al. 2008) at a vacuum wavelength of 4335.2 
A, which corresponds to 4334.0 A in air, a value in close agree- 
ment with our prediction. The significantly stronger /? band of 
hexabenzocoronene occurs near 3344 A (Rouille et al. 2009). 

Unlike the a-band, the /3-band shows features that are unre- 
solved in the matrix spectra. They correspond to numerous bands 
that arise through the interaction of the Sp electronic state with 
the vibrational manifold of S 2 (Rouille et al. 2009). As a con- 
sequence, each peak of the /3-band in the matrix spectra is ex- 
pected to be resolved in the gas phase into a group of bands. For 
example, such groups are found in the spectrum of the S 2 <— So 
transition of pyrene (cf. Sect. 4.1). In the case of pyrene, the ex- 
trapolation from peaks measured by us at 3302 ± 4 and 3234 ± 
5 A in Ar and Ne matrices, respectively, gives an expected posi- 
tion of 3213 ± 8 A of this band in the gas phase. The strongest 
peak of this group has been measured at 3208.2 A or 31 160.9 
+ 0.5 ctrT 1 in the gas phase (Rouille et al. 2004), in good agree- 
ment with the result of the extrapolation. Considering now the 
matrix spectra of hexabenzocoronene, the relative intensities of 
the peaks in the /3-band vary depending on the matrix material. 
We assume that the strongest peak in the Ne matrix spectrum 
corresponds in the gas phase to the group of bands that includes 
the strongest one. The peak lies at 3420 + 4 and 3362 ± 5 A in 
Ar and Ne (Rouille et al. 2009), respectively. Our simple model 
leads to an estimated wavelength of 3344 + 8 A of this band in 
the gas phase. 

3. Observations and data reduction 

We have obtained high-resolution spectra toward the heavily 
reddened supergiants HD 169454 (B1.5Ia, £b-v = 1-13 mag), 
BD -14°5037 (B1.5Ia, £ B -v = 1-56 mag), CPD -32°1734 
(B1.5Ia, £ B _ V = 1-22 mag), CPD -33° 1768 (B2.5I, £ B -v = 
1.46 mag), and HD 183143 (B7Iae, £ B -v = 1-27 mag), using 
the UVES spectrograph at the VLT under programme 079.C- 
0597(A). Total integration times were 14 x 300 sec for CPD 
-32° 1734, 11 x 500 sec for CPD -33° 1768, 16 x 250 sec for 
BD -14°5037, 20 x 50 sec for HD 169454, and 6 x 50 sec for 
HD 183143. The spectra are reduced within the ESO MIDAS 
package following standard procedures and are rebinned to a he- 
liocentric wavelength scale. 

We also use the archival UVES spectra obtained un- 
der the UVES Paranal Observatory Project 'A Library of 
High-Resolution Spectra of Stars across the Hertzsprung- 
Russell Diagram' (Bagnulo et al. 2003) produced under ESO 
Director Discretionary Time (DDT) program 266.D-5655(A). 
The database of Bagnulo et al. (2003) contains a large number 
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Fig. 4. Comparison between the Si «— So origin band of anthracene 
measured in our laboratory (plotted in blue, Staicu et al. (2004) and the 
UVES spectra obtained toward our programme stars. Red lines repre- 
sent synthetic stellar spectra obtained from the models of Gummersbach 
& Kauffer (1993). All stellar spectra have been rebinned to a wave- 
length frame where the detected interstellar absorption lines occur at an 
average radial velocity of km s , compared with the respective rest 
wavelengths in air. 



Fig. 5. Comparison of the weak 5 1 <— 5o transition of phenanthrene (in 
blue) with UVES spectra of our programme stars. See Fig. 4 for further 
details. 



of field stars with spectral types ranging from O to M and lumi- 
nosity classes ranging from main sequence stars to supergiants. 
Here we select a few of the more reddened luminous O- and 
B-type stars with reddenings up to £ B -v = 1 ma g> HD 97253 
(05.5III, £ B -v = 0.50 mag), HD 148937 (06.5, £ B -v = 0.67 
mag), HD 96917 (08.51b, £ B -v = 0.37 mag), HD 76341 (091b, 
£ B -v = 0.46 mag), HD 152003 (09.7Iab, £ B -v = 0.64 mag), 
and toward the reddened B-type supergiants HD 1 12272 (B0.5Ia, 
£ B _v = 0.99 mag), HD 115363 (Blla, £ B _ V = 0.82 mag), HD 
152235 (Blla, £ B _ V = 0.73 mag), and HD 148379 (B 1.51a, 
£ B -v = 0.7 1 mag), Values of £ B _v are from Kazmierczak et al. 
(2010). The archival spectrum of HD 169454 is used for consis- 
tency checks. 

The archive contains several stars with low reddening, which 
we use to identify stellar absorption lines. We select the spec- 
tra toward HD 29138 (Bl lab), HD 58978 (Bl II), and Rigel 
(B8 lab), which are characterized by very faint or absent in- 



terstellar molecular absorption lines of CH 4300 A, CH + 4232 
A, and CN 3875 A, and faint or absent Till 3229 A and Nal 
3302 A absorption lines. We use the models of Gummersbach & 
Kauffer (1993) to assess the stellar continua of our programme 
stars. These models are based on the LTE static plane -parallel 
line-blanketed ATLAS9 model atmosphere of Kurucz (1991) 
and were developed for main-sequence stars with spectral types 
ranging from 09V to B5V. The stars in our sample are super- 
giants in general, some of them with strong stellar winds. Stellar 
parameters such as gravity, rotation, and metallicity which affect 
the width and the depth of the photospheric absorption lines may 
thus deviate significantly from those adopted in the models used 
here. Because of this, the models of Gummersbach & Kauffer 
(1993) may not reproduce the stellar spectra presented here in 
all details, yet the models are still extremely useful to distin- 
guish among stellar and potential interstellar absorption features 
in our UVES spectra. 
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4. Results 

4. 1 . Comparison of astrophysical and laboratory spectra 

In the standard reduction scheme adopted here (cf. Sect. 3), the 
spectra obtained with UVES are rebinned to a heliocentric wave- 
length scale. In order to compare the astronomical spectra with 
our laboratory spectra, we proceed as follows. The astronomical 
spectra contain a number of interstellar absorption lines from di- 
atomic molecules (cf. Sect. 4.3 and Table 4) which we use to 
infer average heliocentric velocities V aV e of -9.8 km s , -7.9 
km s _1 , +36.8 km s _1 , and +35.9 km s _1 of the molecular ma- 
terial toward HD 169454, BD -14°5037, CPD -33°1768, and 
CPD -32° 1734, respectively. We now shift the stellar spectra 
by the respective values of -V ave and plot the modified spec- 
tra, together with the laboratory absorption profiles of the PAHs 
studied here, in Figs. 4-11. Assuming that the sought PAH ab- 
sorbers reside in the molecular material traced by the diatomic 
molecules detected here (cf. Table 4), corresponding interstellar 
absorption features must agree in position and in spectral profile 
with the laboratory profiles plotted in Figs. 4-11. Note that a 
velocity dispersion of up to 4 km s _1 among the various inter- 
stellar absorbers toward a given line of sight would translate into 
a wavelength spread of about one resolution element only, which 
is negligible. 

Figure 4 presents the comparison of the origin band of the 
Si <— So transition of anthracene, which is plotted in blue 
(Staicu et al. 2004), with the UVES spectra of our five pro- 
gramme stars. Red lines represent synthetic stellar profiles ob- 
tained from the models of Gummersbach & Kauffer (1993). 
Toward HD 169454 and BD -14°5037, the expected absorp- 
tion of anthracene agrees in wavelength with a faint photo- 
spheric Felll (3611.7 A) line. Stronger stellar absorption lines 
at 3612.4 A and 3613.6 A correspond to Al III and He I, re- 
spectively. Toward the two B1.5Ia supergiants HD 169454 and 
BD -14°5037, the stellar continuum near the expected absorp- 
tion of anthracene is not very well reproduced by the models of 
Gummersbach & Kauffer (1993), yet the discrepancies are as- 
signed to effects of stellar parameters such as gravity and metal- 
licity. It can thus be concluded that the spectra shown in Fig. 4 
do not reveal signatures of interstellar anthracene absorptions at 
3610.74 A. The observations may nevertheless be used to infer 
upper limits in the column density and in the fractional abun- 
dance of anthracene (cf. Sect. 4.2) toward our lines of sight. 

The laboratory spectrum of the Si <— So origin band of 
phenanthrene at 3409.21 A is given in Fig. 5 together with the 
five UVES spectra. The transition is significantly weaker than 
the S2 <— So origin band, yet the latter occurs at 2826.42 A 
which is not accessible from the ground. The corresponding 
spectral region toward the B7Ia supergiant HD 183143 shown 
in Fig. 5 is dominated by stellar Nill (3407.3 A) and CrII 
(3408.8 A) absorption lines which are well reproduced by the 
stellar models. The hotter atmospheres of HD 169454 and BD 
-14°5037 contain faint absorptions arising from Nell (3406.9 
A) and OH (3407.3 A, 3409.8 A). Toward BD -14°5037, an ab- 
sorption feature occurs near 3408 A which is not reproduced by 
the stellar models. The stellar continuum of the unreddened star 
HD 148379, which is a B1.5 la supergiant as well, is feature- 
less in the spectral region shown in Fig. 5. It is thus tempting to 
assign the absorption feature near = 3408.3 A toward BD 
-14° 5037 to a new diffuse interstellar band, yet its origin from 
phenanthrene is rejected. The absorption feature is significantly 
broader than the phenanthrene band and it would be blue-shifted 
by some 100 km s _1 if due to phenanthrene. BD-14°5037 and 




3344 

wavelength (A) 



Fig. 6. Comparison of the origin band of the Si «— So transition of 
2,3-benzofluorene with the UVES spectra of our programme stars. See 
Fig. 4 for further details. 



HD 169454 are members of the Set OB3 association where high- 
velocity gas has been detected (Federman & Lambert 1992). The 
high- velocity gas, redshifted by some 100 km s _1 , has been de- 
tected toward HD 169454 and a few other stars in the Set OB3 
association, but not toward BD14°5037. Its origin has been as- 
signed to an expanding filamentary shell, probably due to a SN 
explosion in Set OB3. 

Figures 6 and 7 present the comparison of the laboratory 
spectra of 2,3-benzofluorene with our UVES spectra. Figure 6 
displays the vibrationless origin band of the Si <— So transi- 
tion of 2,3-benzofluorene at 3344.16 A. The hotter stars such as 
HD 169454 show faint stellar absorption lines arising from Ne II 
(3344.4 A and 3345.5 A) and Arffl (3344.7 A) which merge 
into a shallow and unresolved absorption feature. The relatively 
featureless stellar continua of our programme stars in this wave- 
length region favor the search of the relatively narrow Si <— So 
transition of 2,3-benzofluorene, yet none of the spectra show 
a corresponding absorption feature near 3344.16 A. Figure 7 
shows a weaker band of 2,3-benzofluorene located at 3266.63 A. 
That band involves an additional vibrational excitation. A strong 
stellar Fe III absorption (3266.9 A) is present in all spectra, and it 
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Fig. 7. Comparison of laboratory data of 2,3-benzofluorene (in blue) 
with UVES spectra of our programme stars. The blue curve represents 
an electronic transition involving the additional excitation of a vibration 
with an energy of 700 cirr 1 (Staicu et al. 2008). See Fig. 4 for further 
details. 



occurs at the expected position of the 2,3-benzofluorene absorp- 
tion. The Fe III line toward HD 169454 exhibits a strong P Cygni 
profile. The adopted models do not reproduce the depth nor the 
profiles of the Felll line toward CPD -32° 1734, HD 169454, 
and BD -14°5037. It can thus not be ruled out that additional 
absorption due to 2,3-benzofluorene occurs near 3266.63 A, yet 
this is unlikely because of the absence of the stronger band at 
3344.16 A. 

Part of the laboratory spectrum of benzo[ghi]perylene show- 
ing two broad vibrational features at 3501.76 A and 3512.15 
A (Tan & Salama 2005) and the corresponding UVES spec- 
tra are shown in Fig. 8. The two vibrational bands belong to 
the 5 2 <— So transition. The coincidence of various stellar ab- 
sorption lines in this wavelength region renders the detection 
of the relatively broad absorption bands of benzo[ghi]perylene 
very difficult. The spectra of the hotter stars show He I absorp- 
tions near 3498.7 A, 3502.4 A, and 3512.5 A, which coincide 
with the absorption wavelengths of benzo[ghi]perylene. Other 
photospheric lines in this region are Nell (3503.6 A) and Cr VII 




3500 3505 3510 3515 

wavelength (A) 



Fig. 8. Comparison of laboratory data of benzo[ghi]perylene (in blue) 
with UVES spectra of our programme stars. The laboratory spectrum 
shows two vibrational bands in the 52 *— So transition (Tan & Salama 
2005). See Fig. 4 for further details. 

(3498.9 A). The cooler stars show Nill (3514.0 A) absorption 
with shallow CrII (3511.8 A) and He I (3512.5 A). The stellar 
models do not reproduce the observed spectra near these wave- 
lengths with sufficient accuracy, and some residuals exist near 
3514 A where benzo[ghi]perylene absorptions are expected. The 
unreddened archival stars of similar spectral type show identi- 
cal photospheric line profiles with broad wings, which leads us 
to reject possible contributions from benzo[ghi]perylene to the 
residuals seen near 3514 A. The strong benzo[ghi]perylene ori- 
gin band of the 52 <— So transition occurs near 3685 A and 
falls into a spectral region which is dominated by strong photo- 
spheric hydrogen absorption lines converging toward the Balmer 
oo — > 2 series limit at 3646 A. Because of the strong photo- 
spheric Balmer lines, a search of the origin band of the S 2 <— So 
transition of benzo[ghi]perylene toward early type supergiants is 
hopeless. 

Our UVES spectra are rather noisy indeed at wavelengths 
below 3210 A where the S 2 <— So absorptions of pyrene occur. 
The signal-to-noise ratio of our spectra is smaller than 100 in 
general, yet the spectra are still useful as the pyrene transitions 
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Fig. 9. UVES spectra in the region where the origin of the 52 <— So 
transition of pyrene occurs (Rouille et al. 2004). See Fig. 4 for further 
details. 



exhibit large oscillator strengths. As mentioned in Sect. 2, the in- 
teraction of the S 2 state of pyrene with the vibrational manifold 
of the S i state gives rise to groups of bands instead of isolated 
features. Figure 9 shows the group of bands of strong intensities 
that appear between 3204 - 3209 A in the place of the origin 
band of S 2 «- S (Rouille et al. 2004). The UVES data reveal 
faint stellar absorption lines due to Felll (3204.8 A) and He I 
(3211.6 A) in the spectra of the hotter stars. The stellar spectra 
reveal no signature that correspond to the group of pyrene bands 
plotted in Fig. 9. Figure 10 compares another group of bands be- 
tween 3163 - 3169 A which correspond to the excitation of an 
in-plane CCC bending mode together with the S 2 electronic state 
of pyrene. These bands coincide with the stellar CII (3165.5, 
3166.0, and 3167.9 A) and Si IV (3165.7 A) absorption lines, 
which makes a detection of pyrene in this wavelength region dif- 
ficult, even in spectra of significantly higher S/N. Thus, the 3204 
- 3209 A region shown in Fig. 9 with its few and faint stellar 
lines appears more promising to search for interstellar pyrene. In 
addition, the pyrene bands in this wavelength region are stronger 
than between 3163 - 3169 A. 




3159 3162 3165 3168 3171 3174 
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Fig. 10. UVES spectra of another region of the S 2 <— So transition 
of pyrene involving electronic and vibrational excitation (Rouille et al. 
2004). See Fig. 4 for further details. 



Figure 11 presents the spectral region where the strongest 
bands of the Sp So transition of hexabenzocoronene is ex- 
pected. The spectral profile for hexabenzocoronene (dashed blue 
line) was measured in a Ne matrix and shifted to the expected 
gas phase position following the method described in Sect. 2.2. 
We are not able to provide absorption cross sections for hexa- 
benzocoronene because our laboratory data does not reveal the 
number of individual gas phase bands that form the broad struc- 
tures seen in Fig. 11. A simulated low-temperature gas phase 
spectrum (solid blue line) that can be compared with the UVES 
spectra is also shown in Fig. 11. It represents two bands at 
the gas phase positions predicted for the strongest peaks in the 
Sp <— So transition of hexabenzocoronene, i.e. at 3344 and 3369 
A. The synthetic bands are represented in terms of Lorentzian 
profiles, where the ratios of their widths and their heights are 
obtained from the spectrum of hexabenzocoronene measured in 
the Ne matrix. We have chosen to fix the widths to 5 and 7 
ctrr 1 because these values are typical of the bandwidths ob- 
served in the jet-cooled spectrum of the S 2 <— So transition 
of benzo[ghi]perylene (Tan & Salama 2005). This transition 
is affected by an interaction between electronic states (Tan & 
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Fig. 11. UVES spectra covering the region where the strongest absorption peaks in the Sp <— So transition of hexabenzocoronene are expected. 
The dashed blue curve represents the hexabenzocoronene spectrum as measured in a Ne matrix and shifted to the expected gas phase position 
following the method described in Sect. 2.2. The two narrow bands, plotted as solid blue curve, are synthetic low-temperature gas phase profiles 
drawn at their predicted gas phase positions (see text). The sharp absorption feature at 3384 A is due to interstellar Till. See Fig. 4 for further 
details. 



Salama 2005) like the Sp <— So transition of hexabenzocoro- 
nene (Rouille et al. 2009). As stated in Sec. 2.2, groups of bands, 
not just a pair, are actually expected in the low-temperature gas 
phase spectrum of hexabenzocoronene in the region covered by 
Fig. 11. Our UVES spectra do not reveal features resembling the 
simulated bands of hexabenzocoronene. 

We have investigated all archival spectra for potential PAH 
absorptions. The archival spectra were obtained toward stars 
with significantly lower reddening, in general, and consequently, 
the S/N ratio in the blue part of the spectrum is higher than in our 
spectra. None of the archival spectra reveal traces of the PAHs 
studied here. This is illustrated in Figs. 12 and 13, which cover 
the spectral region where the two groups of bands of pyrene stud- 
ied here occur. These two pyrene bands have relatively large os- 
cillator strengths (cf. Sect. 2 and Table 1), so relatively small col- 
umn densities of pyrene become detectable. None of the archival 
spectra show signatures of interstellar pyrene, however. The up- 
per limits in the equivalent widths of pyrene are nevertheless 
used to constrain the abundance of pyrene in translucent mate- 
rial, as will be demonstrated in the following section. 



4.2. Constraints of PAH column densities in translucent 
material 



Our UVES spectra together with the laboratory data may 
be used to infer upper limits in the column densities 
of anthracene, phenanthrene, pyrene, 2,3-benzofluorene, and 
benzo[ghi]perylene. In cases where a sought absorption band is 
not detected, an upper limit in the equivalent width Wa of the 
absorption band is given by W A = 3 VMAA/(S/N), where AA 
is the spectral width of a resolution element, M is the number 
of resolution elements over which the sought absorption band 
spreads, and S/N is the signal-to-noise ratio near the expected 
absorption feature (see Bohlin et al. (1983) for a comprehensive 
treatment of error propagation in absorption line measurements). 
The equivalent width W,j of an interstellar absorption line is de- 
fined in terms of W,\ = J (F c - Fj)/F c dA where F A is the mea- 
sured flux level at wavelength A and where F c is the flux in the 
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continuum adjacent to the absorption line. The equivalent width 
is related to the optical depth t(A) according to 

Wi = J (1 - e- TW )dA, (4) 

which, for low optical depths t(A) « 1, is expressed as 

W A = j T(A)dA = N j cr(A)dA. (5) 

N is the number of absorbers per cm 2 and cr(A) is the absorption 
cross section. 

Our UVES spectra were obtained with a spectral resolution 
of R = A/AA = 80000. For simplicity, we adopt a constant 
width of AA = 45 mA per resolution element over the entire 
wavelength range studied here. Table 3 summarizes the central 
wavelengths and the corresponding values of 3 yfMAA for each 
of the sought absorption bands, and the corresponding signal-to- 
noise ratio S/N, the upper limits in Wa, and the upper limits in 
the column densities N toward HD 169454, BD -14°5037, CPD 
-33°1768, CPD -32°1734, and HD 183143. 

Our data indicates maximal column densities of anthracene, 
pyrene, and benzo[ghi]perylene of a few 10 12 cirT 2 toward the 
heavily reddened programme stars of our UVES sample. Toward 
the stars retrieved from the archive, upper column densities in 
pyrene are below values of 10 12 crrT 2 . Using the relations of 
Mattila (1986), the upper limits in the PAH column densities 
and the observed CH column densities (cf. Sect. 4.3) are used 
to infer upper limits in the fractional abundances of /(PAH) = 
AT(PAH)/7Y(H) = 1.1 x 10- 8 7Y(PAH)/7Y(CH). Our observations 
indicate values of /(anthracene) < 1.5 - 2 x 10" 10 toward 
HD169454, BD-14°5037, CPD-32°1734, and CPD-33°1768, 
and /(anthracene) < 10" 9 toward HD 183143. The upper lim- 
its in the fractional abundances of pyrene and 2,3-benzofluorene 
are < 5 x 10 10 toward HD 169454 and BD-14°5037 and 
< 3.5 x 10~ 9 toward the remaining three stars. Fractional abun- 
dances of pyrene toward the more diffuse lines of sight traced by 
the archival stars are below a few times 10~ 10 . For phenanthrene, 
which has a relatively small oscillator strength in the transitions 
studied here (cf. Table 1), upper limits in the fractional abun- 
dances are a few 10~ 8 . Our values may be compared to fractional 
abundances relative to H2 derived for other complex molecules 
from radio telescope observations. Marcelino et al. (2007) report 
a relative abundance of 4 x 10~ 9 for propylene detected toward 
the dark cloud TMC-1 while Belloche et al. (2009) derived abun- 
dances for ethyl formate and n-propyl cyanide observed toward 
Sgr B2(N) of 3.6 x 10~ 9 and 1.0 x 10~ 9 , respectively. 

4.3. Interstellar molecular absorption lines 

Our UVES spectra contain various interstellar absorption lines 
which are presented for consistency checks with previous re- 
sults. 

OH + . Our spectrum obtained toward CPD -32° 1734 con- 
tains an absorption line near 3584 A which we assign to the in- 
terstellar R n (0) (i air = 3583.769 A) line of the (0,0) band of 
the A 3 H - X 3 Z~ system of OH + . The relevant part of the spec- 
trum is reproduced in Fig. 14, together with the spectra of HD 
169454 and BD -14°5037. Using the band oscillator strength of 
/oo = 3.0764 x 10~ 3 given by de Almeida & Singh (1981), we 
convert the measured equivalent width of W A =4 mA into an 
OH + column density of N(OU + ) = 11 x 10 12 cm 2 . The OH + 
absorption line appears at a heliocentric velocity of + 35.4 km 



s 1 which is consistent with the average velocities of interstellar 
CN (V hel = +37.5 km s" 1 ), CH + (V hel = +35.6 km s" 1 ), and CH 
(V hel = +35.6 km s" 1 ) toward that line of sight (cf. Table 4). 

The OH + radical has only recently been discovered in the 
interstellar medium. Wyrowski et al. (2010) detected the OH + 
N - 1-0, 7 = 0-1 absorption line toward the strong con- 
tinuum source Sgr B2 using APEX at sub-mm wavelengths, and 
Gerin et al. (2010) report the detection of OH + absorption to- 
ward G10.6-04 (W31C) using HIFI on HERSCHEL. Optical ab- 
sorption lines arising from interstellar OH + have been searched 
by de Almeida & Singh (1981) without success. Optical detec- 
tions of OH + are reported by Krelowski et al. (2010a). Using the 
relations of Mattila (1986), the observed columns of OH + and 
CH are used to infer a fractional abundance of N(OH + )/N(H) = 
1.1 x 10- 8 JV(OH + )/jV(CH) = 9 x 10- 10 , where N(H) = 2N(H 2 ). 
The fractional abundance of OH + of about 10~ 9 cirT 2 is almost 
two orders of magnitude lower than the abundances inferred to- 
ward the massive star forming region W31C (Gerin et al. 2010). 
The absence of velocity shifts between OH + and the neutral in- 
terstellar molecules (cf. Table 4) argues against a production of 
OH + in shocks, as proposed by de Almeida (1990). Interstellar 
shocks have been ruled out as the general production site of CH + 
(Gredel et al. 1993). 

NH. Interstellar NH has been discovered by Meyer & 
Roth (1991). From their measurements of the R\(Q) (A-^ = 
3358.0525 A) and R Q 21 (0) (A. dir = 3353.9235 A) lines of the 
(0,0) band of the A 3 n - X 3 E system, the authors inferred col- 
umn densities of 1 x 10 12 cirT 2 and 2.7 x 10 12 cirT 2 toward £ 
Per and HD 27778, respectively. Both lines are detected in our 
spectrum of HD 169454 and are marginally present toward CPD 
-32° 1734. We use oscillator strengths of / = 0.0041 for the 
Ri(0) and / = 0.0024 for the R Q 2 i(0) line (Meyer & Roth 1991) 
to infer a column density of iV(NH) = 7.6 x 10 12 cirT 2 toward 
HD 169454. The column density is in good agreement with the 
value of 6.7 X 10 12 cm" 2 given recently by Weselak et al. (2009a) 
in their comprehensive analysis of interstellar NH. An NH col- 
umn density of JV(NH) « 2.4x 10 12 cirT 2 is inferred toward CPD 
-32° 1734. Toward both lines of sight, the radial velocities of the 
NH absorption lines agree with those of CH, CH + , and CN. 

CN. The R(l) (A. dil = 3579.453 A), R(0) (A. dil = 3579.963 A), 
and P(l) (i air = 3580.937 A) lines of the (1,0) band of the 
B 2 S - X 2 S system of CN are detected toward HD 169454, and 
the R(0) line is detected toward BD -14°5037, CPD -32° 1734, 
and CPD 33° 1768. We use the molecular parameters given by 
Meyer et al. (1989) to infer interstellar CN abundances of N( J = 
0) = 2.4xl0 13 cirr 2 andjV(/= 1)= 1.2 x 10 13 cm" 2 toward HD 
169454, in the limit of Doppler-Z? values of b — > 00. The ratio in 
the population densities corresponds to an excitation tempera- 
ture of r ex = 2.9 K. The value of N(J = 0) inferred here is con- 
sistent with the value of N(J = 0) = 3 x 10 13 cirT 2 obtained for a 
Dopper-fe value of b vr = 0.45 km s _1 (Gredel et al. 1991), where 
b vr was obtained from an analysis of CN absorption lines arising 
in the violet and red systems. The CN column densities in J = 
toward BD -14°5037, CPD -32° 1734, and CPD -33° 1768 (cf. 
Table 4) are consistent with the measurements of Gredel et al. 
(1991) and Gredel et al. (2002) toward those lines of sight. 

CH. We report the detection of the R 2 (l) (i air = 3137.576 A) 
and p Qn(l) i^k = 3145.996 A) lines and of the [Q 2 (l) + Q 
Rn(l)] (/lair = 3143.183 A) unresolved line blend of the (0,0) 
band of the C 2 Z + -X 2 n system of CH. In addition, the [ G /?i 2 (l)+ 
Gi 2 (l)] (^air = 3633.29 A) line blend and the p Q l2 {\) (A mr = 
3636.27 A) line of the (1,0) band of the B 2 E~ - X 2 U system 
are detected. The measured equivalent widths of the lines aris- 
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Table 3. Upper limits of large molecules toward the lines of sight studied here. 





HD 169454 




BD - 


14°5037 




CPD 


-33° 1768 


CPD 


-32° 1734 


HD 183143 




3 y/MAA 


S/N 


w A 


N a 

max 


S/N 


w A 


N" 

max 


S/N 


W A 


N" 

max 


S/N 


W A 


N" 

' max 


S/N 


W A 


N a 

1 y max 


anthracene 
































3610.74 371 


230 


1.6 


0.8 


113 


3.3 


1.7 


81 


4.6 


2.4 


143 


2.6 


1.4 


69 


5.4 


2.8 


phenanthrene 
































3409.21 284 


165 


1.7 


105 


97 


2.9 


180 


62 


4.6 


285 


145 


1.9 


118 


53 


5.3 


330 


pyrene 
































3208.22 1493 


70 


21 


2.3 


37 


40 


4.4 


23 


65 


7.2 


51 


29 


3.2 


23 


65 


7.1 


2,3-benzofluorene 
































3344.16 270 


117 


2.3 


2.7 


88 


3.1 


3.7 


47 


5.7 


6.8 


117 


2.3 


2.8 


41 


6.6 


7.9 


benzo [ghi]perylene 
































3512.15 1822 


206 


8.8 


8.8 


125 


14.6 


14.6 


82 


22 


22 


165 


11 


11 


76 


24 


24 


3501.76 1465 


206 


7.1 


12.8 


125 


11.7 


21 


82 


18 


32 


165 


8.9 


16 


76 


19 


34 



" Upper limits N miix in units of 10 12 cm 2 



ing in the C-X system are converted into column densities using 
the molecular parameters and the method of Lien (1984). For 
the oscillator strength of the (1,0) band of the B-X system, we 
scale the value of the (0,0) band of Lien (1984) using the ratio of 
7.8 of the band oscillator strengths of the (0,0) and (1,0) bands 
as inferred from the lifetime measurements of Brooks & Smith 
(1974). The observed CH absorption lines and line blends arise 
from the lower or upper component of the J = 1 /2 A-doublet. 
In Table 4, we list total CH column densities N to t, which are ob- 
tained under the assumption that N tot = 2N\, 2N U , where Nj and 
N u are the column densities in the lower and upper A compo- 
nents, respectively. For the purpose of the present consistency 
check, we assume that all lines are optically thin, and ignore a 
detailed curve of growth analysis. 

CH + . The R(0) line (i air = 3745.307 A) of the (2,0) band of 
the A'n - X : 2 + system of CH + is detected toward all five lines 
of sight studied here, and the R(0) line of the (3,0) band (/l air = 
3579.021 A) is detected toward HD 169454, BD -14°5037, 
CPD -32° 1734, and CPD -33° 1768. We use the molecular 
parameters given by Weselak et al. (2009b) to infer CH + col- 
umn densities ranging from N(CH + ) = 1.6 X 10 13 crrT 2 toward 
HD 169454 to N(CR + ) « 4 x 10 13 ctrT 2 toward BD-14°5037, 
CPD-32°1734, and CPD-33°1768. 
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5. Discussion 

5.1. The long and eventful history of DIB candidates 

Over the last decades, a large variety of potential DIB carriers 
has been proposed, either because laboratory spectra revealed 
coincidences with absorption features in astronomical spectra or 
because theoretical calculations predicted molecular transition 
wavelengths that matched those of observed DIBs. For instance, 
Tulej et al. (1998) proposed based on a few coincidences 
between bound-bound electronic transitions of and five DIB 
absorption wavelengths. The proposition is no longer generally 
accepted (McCall et al. 2001), as the exact wavelengths and ab- 
sorption profiles of the lines do not match the interstellar 
profiles closely enough. A number of close agreements in the 
absorption wavelengths of neutral carbon chain molecules and 
their cations with DIB absorptions were found by Motylewski 
et al. (2000), yet in all cases small but significant shifts be- 
tween the laboratory and the interstellar absorption features re- 
mained. These shifts could not be explained by radial veloc- 
ity shifts of the interstellar material nor could they be traced 
to calibration issues in the astronomical data. The HC4H + di- 
acetylene cation has recently been purported as a DIB carrier 
(Krelowski et al. 2010b). The authors averaged a large number 
of stellar spectra of different stars in order to improve the S/N- 
ratio. The average spectrum reveals a faint absorption feature 
near 5069 A which Krelowski et al. (2010b) assign to a previ- 
ously undetected diffuse interstellar band. The gas-phase spectra 
obtained by Motylewski et al. (2000) shows a prominent absorp- 
tion band of HC/tH + at 5069 A, which matches exactly the po- 
sition of the 5069 A feature. However, the HC4FF band profile 
obtained from the laboratory data deviates significantly from the 
5069 A feature, and it is difficult to see how a convolved HC/jH + 
band profile, assuming typical Doppler-Z? broadening of the ro- 
tational lines and macroscopic turbulent gas motions in translu- 
cent clouds, can be brought into agreement with the profile of the 
5069 A absorption feature reported by Krelowski et al. (2010b). 

In a similar fashion, the recent evidence provided by Maier 
et al. (2010) that l-C 3 U 2 is responsible for the 5450 A and 4881 
A DIBs requires some scrutiny. While the authors convincingly 
assign two absorption features near 4887 A and 5450 A seen in 
their laboratory CRDS spectra to the 2 l Q and 2^ transitions in the 
B'Bj <— X X A\ system of Z-C3H2, the agreement with interstellar 
DIBs is good for the 5450 A line only. The 4887 A transition of 
Z-C3H2 does not match the exact shape and position of the 4881 
A DIB, unless it is assumed that several DIB carriers contribute 
to this feature. More critically, Ne-matrix spectra also reveal an 
Z-C3H2 absorption feature which arises from the 2 ( ,4g transition, 

which is about a factor of 2 - 3 fainter than the 4887 A and 5450 
A transitions of /-C3H2. The strength of the 2^4,, transition is be- 
low the detection limit of the CRDS experiments presented by 
Maier et al. (2010) and an accurate gas phase wavelength is not 
available, yet its position is predicted to be near 5173 A. Maier 
et al. (2010) did search their spectra for an interstellar absorp- 
tion band near 5173 A, yet they note that the presence of vari- 
ous stellar absorption lines hampered their effort. The high S/N 
spectrum presented by Hobbs et al. (2009) toward HD 183143 
does show broad and prominent DIB absorptions near 4881 A 
and 5450 A. The latter spectrum is relatively clean in the 5165 - 
5180 A region, yet it does not show the expected signature from 
the 2 ( '4g transition of /-C3H2 in this wavelength interval. 

Other molecular cations and anions have served as prominent 
candidates of DIB carriers in the past. Matrix isolation experi- 



ments of Cg revealed absorption features near the 9577 A and 
9632 A DIBs (Foing & Ehrenfreund 1997), however a confirma- 
tion of the Cg gas phase wavelengths based on CRDS spectra 
is not yet available. That work supported the general proposition 
that PAH cations are responsible for DIB absorptions (Salama 
et al. 1996, 1999). Jet-cooled absorption spectra measured in the 
laboratory (Romanini et al. 1999; Biennier et al. 2003, 2004; 
Sukhorukov et al. 2004; Tan & Salama 2006) were subsequently 
obtained, and interstellar spectra obtained toward the Perseus 
molecular cloud complex did reveal weak features that coincide 
with absorption bands of the naphthalene and anthracene cations 
(Iglesias-Groth et al. 2008, 2010). In the case of the anthracene 
cation, the CRDS spectrum of Sukhorukov et al. (2004) reveals 
an absorption at 7085.7 A which has a FWHM of 47 A. While 
the spectrum of Iglesias-Groth et al. (2010) may show an inter- 
stellar absorption band near that wavelength, the width of the 
purported interstellar absorption feature is significantly smaller 
than 47 A. 

Potential DIB carriers have also been proposed from the- 
ory. With rotational contour fitting, Kerr et al. (1996) were 
able to reproduce the structure seen in the 6614 A and 5797 A 
DIBs assuming that they are due to planar oblate symmetric top 
molecules, such as large carbon ring molecules with 14 - 30 
atoms, rather than PAHs. An excellent agreement was found be- 
tween the calculated absorption wavelength of the 0[j band of the 

l B x - X l A' electronic system of CH 2 CN and the 8037 A DIB 
(Cordiner & Sarre 2007). However, the calculated absorption 
wavelength needs to be confirmed by gas phase laboratory ex- 
periments. It is also noted that the match between the calculated 
CH2CN absorption and the 8037 A DIB is good for an assumed 
rotational excitation temperature of 2.7 K and a Doppler-Z? value 
of 16 - 33 km s _1 . At smaller ^-values, the rotational structure of 
CH 2 CN should be resolved in the astronomical spectra, in dis- 
agreement with the observed broad and diffuse absorption near 
8037 A. It is also noted that at a nominal ortho to para ratio of 
3:1, the calculated spectrum of CH2CN does not match the ob- 
served spectra. Similarly, a recent claim that proflavine and its 
anion, cation, and di -cation are responsible for the 4066 A, 4363 
A, 4175 A, and 5259 A DIBs (Bonaca & Bilalbegovic 2010) re- 
quires verification. The proposition is based on computed opti- 
cal absorption wavelengths using pseudo-potential density func- 
tional theory methods. The matches in wavelengths are marginal 
at most, and need to be confirmed by laboratory experiments. 
In any case, as pointed out at the end of Sect. 2.1, extreme care 
must be taken if computational results on electronic transitions 
are used for the interpretation of astronomical observations, due 
to the large uncertainties still being involved in present-day cal- 
culations. 



5.2. Fractional PAH abundances in translucent molecular 
clouds 

PAHs are believed to be formed in the outflows of carbon-rich 
AGB stars from where they are deposited into the interstellar 
medium. The fractional PAH abundance is estimated to be of the 
order of /(PAH) = N(PAU)/N(U) = 10~ 7 (Habart et al. 2004), as- 
suming some 100 carbon atoms per PAH molecule. In the diffuse 
interstellar material, at densities of «h = 100 cm" 3 , a temperature 
of T = 100 K, and fractional electron abundances of n e = 10" 4 , 
small PAHs such as anthracene, phenanthrene, and pyrene, are 
rapidly destroyed by the ambient UV radiation (Le Page et al. 
2003). The photoionization of PAHs in diffuse clouds is a ma- 
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jor heating source (Lepp & Dalgarno 1988a) and supersedes the 
grain heating for PAH abundances greater than 2 x 10~ 7 . As 
the visual extinction increases, the abundance of PAHs increases 
rapidly, as shown by Lepp et al. (1988). The authors postulated 
a fractional PAH abundance of 10~ 7 in an attempt to resolve 
density discrepancies which persisted for translucent molecular 
cloud models. The transition between neutral PAHs and their an- 
ions and cations is governed by the balance between photoion- 
ization and recombination as well as electron attachment and 
photo-detachment. At visual extinctions of Ay = 1 mag (which 
corresponds to a reddening of £b-v = 0.3 mag), the fractional 
abundances of neutral PAH and their anions are comparable and 
exceed the fractional abundance of the positively charged PAHs 
by about two orders of magnitude. As the visual extinction in- 
creases further, PAHs form a drain to the free electrons for PAH 
abundances greater than 1CT 8 , resulting in a substantial reduc- 
tion in the density of electrons (Lepp & Dalgarno 1988b). As 
a result, the molecular ions undergo nondestructive neutraliza- 
tion reactions with the negatively charged PAHs, which leads to 
a significant increase of the fractional abundances of the carbon- 
bearing molecules. 

Based on the infrared emission bands, Allamandola et al. 
(1985) evaluated the fractional abundance of interstellar PAHs 
to be 2 x 1(T 7 . Later, Allamandola et al. (1989) distinguished 
between small and large PAHs, which were defined as contain- 
ing 25 and 300-400 C atoms, respectively. A minimal fractional 
abundance of 3 x 10~ 8 was obtained for the small PAHs while a 
value of 8 x 10~ 9 was derived for the large PAHs (Allamandola 
et al. 1989), corresponding to fractions of a 1 - 5% and 0.5% 
of elemental carbon locked in small and large PAHs, respec- 
tively. Considering a single class of PAHs with a typical size 
of 50 C atoms, an abundance of about 3 x 10~ 7 was derived by 
Tielens (2008). Of interest to our study is the minimal interstellar 
fractional abundance of the small PAHs. With the exception of 
hexabenzocoronene, the molecules we have studied can be con- 
sidered to belong to the class of small PAHs. The upper limits we 
have derived for the fractional abundances of anthracene, pyrene, 
and 2,3-benzofiuorene are about two to three orders of magni- 
tude lower than the estimates presented by Tielens (2008). The 
fractional abundance of phenanthrene is not well constrained 
from our observations, mainly due to the fact that the phenan- 
threne band studied here has a small oscillator strength. We 
would thus conclude that small and neutral PAHs must occur in 
a large variety (several 100 species) in translucent clouds, if it is 
assumed a fraction of several percent of the elemental carbon is 
locked up in small PAHs. Alternatively, small and neutral PAHs 
may not be very abundant in translucent material, either because 
small PAHs carry positive or negative charge or because signifi- 
cantly less elemental carbon is locked up in these molecules. 

In their investigation of the ultraviolet interstellar extinction 
curve using the HST/STIS spectrograph, Clayton et al. (2003) 
did not detect absorption bands in the 1150 A - 3180 A wave- 
length region that could be associated with PAHs. The authors 
conlcuded that either PAHs are destroyed by the interstellar ra- 
diation field, or that they occur in such a large variety that indi- 
vidual absorptions are too weak to be detected. 

Concerning the larger PAHs, an upper limit in the column 
density of hexabenzocoronene of 4 x 10 12 cm" 2 in translucent 
molecular material has recently been estimated by Kokkin et al. 
(2008) from spectra covering the vibronic hexabenzocoronene 
absorptions in the a-band near 4262 A. The a-band is, how- 
ever, not well suited for the search of interstellar hexabenzocoro- 
nene, as it is significantly weaker than the /3-band near 3344 A 



(Rouille et al. 2009). Unfortunately, we are not able to constrain 
the molecular abundance of hexabenzocoronene from our obser- 
vations, as we were unable to derive absorption cross sections of 
the jfi-band (cf. Sect. 2.1 and Sect. 4). We nevertheless propose 
that searches for interstellar hexabenzocoronene be conducted 
near 3344 A. 



5.3. Molecular anions in translucent material 

Carbon-chain molecular anions have recently been established 
as a major constituent of dense material. The detection of CeH 
in the TMC-1 molecular cloud by McCarthy et al. (2006) marks 
the discovery of molecular anions in the interstellar medium. 
That work was followed by the observation of several other 
anions in interstellar and circumstellar material, e. g. C4H 
(Cernicharo et al. 2007; Agundez et al. 2008), C 8 H (Remijan 
et al. 2007), C 3 N~ (Thaddeus et al. 2008), and CN (Agundez 
et al. 2010). The effects of molecular anions on the chemistry of 
dark clouds have been rediscussed by Millar et al. (2007) and 
Walsh et al. (2009). The models set up for dark clouds produce 
large amounts of molecular anions at molecular cloud ages be- 
low 10 5 yr. The anion abundances are proportional to the elec- 
tron abundance and are characterized by a sharp drop by several 
orders of magnitude as the clouds reach ages of 10 6 yr or above. 
The anion to neutral ratio is typically of the order of a few per- 
cent. 



5.4. Individual lines of sight 

HD 169454 and BD -14°5037. Various studies of the lines of 
sight toward HD 169454 and BD -14°5037 are available in the 
literature. The C2 observations of Gredel & Muench (1986) re- 
vealed the presence of very cold gas toward HD 169454, with 
gaskinetic temperatures of 7\i n = 15 K and gas densities of n = 
400 cm 4 . The C2 absorption occurs at heliocentric velocities of 
V hel = -8.5 km s _1 . Toward BD -14°5037, the C 2 absorption oc- 
curs in two velocity components at Vh e i = -6.0 km s _1 and Vh e i 
= + 4.0 km s -1 , where the gas in the two velocity components is 
characterized by temperatures and densities of 7\i n = 30 K and 
n = 400 cirT 3 and T^, = 50 K and n = 250 cm 4 , respectively. 
The CH and CN absorption lines are resolved into two velocity 
components at similar heliocentric velocities (cf. Table 4). The 
CH + absorption profiles are unresolved and appear near helio- 
centric velocities which are consistent with those of C2, CH, and 
CN. Very small Doppler-fo values of b = 0.4 km s _1 toward both 
lines of sight were inferred by Gredel et al. (1991) from a com- 
parison of CN column densities derived from the (0,0) bands of 
the CN red and violet systems. A curve of growth analysis us- 
ing b = 0.4 km s _1 resulted in total CN column densities of 
AT(CN) = 4.2 x 10 13 cm 4 and AT(CN) = 1.6 X 10 13 cirT 2 to- 
ward HD 169454 and BD-14°5037, respectively, and in very 
good agreement with the present results. The present data are 
obtained from absorption lines in the (1,0) band, which are about 
a factor of 10 weaker than the absorption lines in the (0,0) band. 
The excellent agreement between both sets of measurements 
confirms the very small Doppler-Z? values toward both lines of 
sight. Using the small Doppler-Z? values, Gredel et al. (1993) ob- 
tained CH and CH + column densities of JV(CH) = 4.4 x 10 13 
cm 4 and N(CH + ) = 1.8 x 10 13 cirr 2 toward HD 169454, and 
JV(CH) = 7.3 x 10 13 cm" 2 and N(CH + ) = 6.4 x 10 13 cnr 2 toward 
BD -14°5037, again in excellent agreement with the present re- 
sults. 
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CPD -32° 1734 and CPD -33° 1768. Very large column den- 
sities of C2, CH, CH + , and CN have previously been measured 
toward the heavily reddened supergiants CPD -32° 1734 and 
CPD -33° 1768 in the NGC 2439 association (Gredel et al. 2002, 
and references therein). From observations of C2, gaskinetic 
temperatures of 7\i n = 85 K, densities of n > 1000 cm" 3 , and 
column densities of N(C2) ~ 10 14 cm" 2 were inferred toward 
both lines of sight. The CH and CH + column densities reach 
values of JV(CH + ) = 4.4 x 10 13 cm" 2 and /V(CH) = 9 x 10 13 
cm" 2 toward CPD -32° 1734 and N(CH + ) = 6.0 x 10 13 cm" 2 and 
N(CH) = 12 X 10 13 cm" 2 toward CPD -33°1768. The column 
densities and heliocentric velocities obtained here are in excel- 
lent agreement with the previously published results. 

HD 183143. A comprehensive optical study of the diffuse in- 
terstellar bands toward HD 183143 has recently been introduced 
by Hobbs et al. (2009). The authors present a catalog of 414 
DIBs from spectra obtained in the 3900 - 8100 A wavelength 
region with a signal-to-noise ratio of S/N as 1000 and a spectral 
resolution of R — 38 000. Their study forms part of a new sur- 
vey of DIBs toward 30 OB stars carried out at the Apache Point 
Observatory at a very high signal-to-noise ratio. Molecular car- 
bon was not detected toward HD 183143 (Gredel 1999), which 
indicates that the line of sight is characterized by low-density 
material. 

6. Summary 

We have presented high signal-to-noise (S/N > 100) absorption 
line spectra toward heavily reddened early type supergiants in 
the 3050 - 3850 A wavelength region. The spectra are compared 
with laboratory gas-phase spectra of anthracene, phenanthrene, 
pyrene, 2,3-benzofluorene, and benzo[ghi]perylene. We have de- 
veloped methods to infer absolute absorption cross sections from 
the laboratory spectra, which we use to derive upper limits in 
the column densities and the abundances of these neutral PAHs. 
Our methods lead to reasonable values for the absorption cross 
sections, if compared with cross sections derived from ab initio 
methods. The fractional abundances for the PAHs are below a 
few times 10" 10 for anthracene, pyrene, and 2,3-benzofluorene. 
Fractional abundances of benzo[ghi]perylene below a few times 
10" 9 and of phenanthrene below a few times 10" 8 are derived. 
Thus, a large variety of small and neutral PAHs such as an- 
thracene and pyrene is required to explain a canonical value of 
10" 7 for the total fractional PAH abundance in the interstellar 
medium. Alternatively, it may be concluded that small and neu- 
tral PAHs are suppressed in translucent molecular clouds, either 
because the small PAHs carry predominantly positive or nega- 
tive charge or because they don't survive the prevailing physical 
conditions in such clouds. 
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